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ABSTRACT 

Transmission spectroscopy at UV wavelengths is a rich and largely unexplored source of information about 
the upper atmospheres of extrasolar planets. So far, UV transit observations have led to the detection of atomic 
hydrogen, oxygen and ionized carbon in the upper atmosphere of HD209458b. The interpretation of these 
observations is controversial - it is not clear if the absorption arises from an escaping atmosphere interacting 
with the stellar radiation and stellar wind, or the thermosphere inside the Roche lobe. In this paper, we introduce 
an empirical model that can be used to analyze UV transit depths of extrasolar planets. We use this model to 
interpret the transits of HD209458b in the H Lyman a and the 1304 A O I triplet emission lines. The results 
indicate that the mean temperature of the thermosphere is J = 8,000-1 1,000 K and that the H2/H dissociation 
front is located at pressures between p = 0.1-1 /ibar, which correspond to an altitude of z ss 1.1 R^,. The upper 
boundary of the model thermosphere is located at altitudes between z = 2.7-3 R^, above which the atmosphere 
is mostly ionized. We find that the H I transit depth reflects the optical depth of the thermosphere in the wings 
of the H Lyman a line but that the atmosphere also overflows the Roche lobe. By assuming a solar mixing ratio 
of oxygen, we obtain an O I transit depth that is statistically consistent with the observations. An O I transit 
depth comparable to or slightly larger than the H I transit depth is possible if the atmosphere is undergoing fast 
hydrodynamic escape, the O/H ratio is supersolar, or if a significant quantity of neutral oxygen is found outside 
the Roche lobe. We find that the observations can be explained solely by absorption in the upper atmosphere 
and extended clouds of ENAs or atoms strongly perturbed by radiation pressure are not required. Due to the 
large uncertainty in the data, repeated observations are necessary to better constrain the O I transit depths and 
thus the composition of the thermosphere. 

Subject headings: ultraviolet: general — plasmas — hydrodynamics — planets and satellites:general — 
ISM:lines and bands 



1. INTRODUCTION 

Visible and near-IR spectroscopy of transiting extrasolar 
planets has led to the dete ction of Na I, H2O, CH4 , CO, and 
CO2 (Charbonne au et al . 2002; Tin etti et al.][2007l [Swain et| 
ir. 2008 , 2009 ) in the atmospheres of extrasolar giant planets 
tEGPs). Discoveries such as these have helped to shift the em- 
phasis in the study of extrasolar planets from mere detection 
towards characterization of the dynamics and composition of 
their atmospheres. Modeling and observations of a large va- 
riety of objects allow us to constrain theories of the forma- 
tion and evolution of different planets and their atmospheres. 
In this respect, studies of mass loss from the atmospheres of 
close-in EGPs and terrestrial planets are of particular inter- 
est. In addition to being interesting in their own right, these 
studies are valuable for constraining models of planetary at- 
mospheres during the early history of the solar system. 

Approximately 25% of the currently known extrasolar plan- 
ets orbit within 0. 1 AU from their host star^ At such small 
orbital distances many of these planets are subject to intense 
stellar in^adiation. In particular, the outermost layers of the 
atmosphere are strongly ionized and heated by the absorp- 
tion of FUV and EUV radiation, which leads to the formation 
of the thermosphere ( e.g., Yelle 2004 ; Garcia Munoz 2007j 
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spheres of close-m extrasolar giant planets are hot, greatly 
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extended, and their composition is dominated by atomic con- 
stituents. The same should apply to close-in super-Earths or- 
biting M dwarf stars. These are the closest types of object 
to an Earth-like habitable planet whose atmospheres we may 
be able to probe with space-b orn or even ground-base d tele- 
scopes in th e near future (e.g., |Palle et al . 2009; Charbo nneau| 
et al. [200^ 1. As a consequence, the UV transit depths of ex- 
trasolar planets should be much larger than the visible and IR 
transit depths because the observations probe the upper atmo- 
sphere that strongly absorbs UV radiation in electronic reso- 
nance lines. 

Transit observations of extrasola r plane ts in the FUV were 
pioneered by [Vidal-Mad jar et al J ( [2003) (hereafter VM03) 
and [Wdal-Madj ar et al.|Ep04| l (hereaFter VM04) who de- 
tecteaatomic hydrogen (H I), oxygen (O I), and ionized car- 
bon (C II) in the upper atmosphere of the transiting planet 
HD209458b. All of these detections were obtained by using 
the Space Telescope Imaging Spectrograph (STIS) onboard 
the Hubble Space Telescope (HST). VM03 used the G140M 
grating to observe three transits of HD209458b in the wave- 
length range covering the stellar H Lyman a (H Lya) emis- 
sion line. The spectral resolution of these observations was 
^0.08 A, which allowed for the details of the line profile to 
be resolved. VM04 used the G140L grating with a low spec- 
tral resolution of ^2.5 A to observe four transits in the wave- 
length range of [1180,1710] A and detected absorption in the 
O I and C II stellar emission lines. 

VM03 deduced a H I transit depth of 15 ± 4 % from 
the ratio of the flux in two wavelength regions around the 
core of the H Lya line to the flux in the wings of the line 
during transit. Based on this observation, they concluded 
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that H I overflows the Roche lobe with a mass loss rate of 
dM /dt > 10^ kg s"'. They also suggested that the planet is 
followed by a cometary tail that is shaped by stellar radia- 
tion p ressure acting on the escapin g hydrogen . Late r, |Ben-| 
|J^r ("2007 ) (hereafter BJ07) and Ben-Jaffel (lOOS") (hefe^ 
after BJ08) presented a thorough and convincing reanalysis of 
the same data. Based on this analysis, BJ08 argued that there 
is no evidence for a cometary tail in the transit light curve or 
the detailed transit depth measurements. Further, he showed 
that the observed absorption may arise from atomic hydrogen 
below the Roche lobe and that the H I absorption line profile 
is broadened by thermal and natural broadening in the ther- 
mosphere of the planet. However, his results also imply that 
hydrogen overflows the Roche lobe, and thus the atmosphere 



is still evaporating (IVidal-Madjar et al.|2008[ ). 

VM04 used the low resolution G140L data to obtain O I 
and C II transit depths of 13 ± 4.5 and 7.5 ± 3.5 %, respec- 
tively. They also calculated a full-width H I transit depth of 
5 ± 2 % from the unresolved H Lya line and claimed that 
this transit depth is consistent with the much stronger absorp- 
tion observed within a limited section of the line profile by 
VM03. Further, they argued that the large O I and C II transit 
depths are possible because these species overflow the Roche 
lobe and the absorption lines are broadened by the velocity 
dispersion of the escaping gas. Thus the observations were 
interpreted as proof that the atmosphere is undergoing fast hy- 
drodynami c escape. 

Ben-Jaff eTand Hosseirii| ( |2010| l (hereafter BJIO) pubHshed 
a thorough reanalysis of the G140L data. They obtained 
revised full-width H I, O I, and C II transit depths of 
6.6 ± 2.3 %, 10.5 ± 4.4 %, and 7.4 ± 4.7 %, respectively. 
We note that these depths are only 2, 1.93, and 1.15(t, respec- 
tively, away from the FUV continuum transit depth of ^2 %. 
The more detailed H I transit depth measurements reported 
by BJ08 provide stronger constraints for the H Lya line but 
no such constraints are available for the O I and C 11 tran- 
sits. In order to explain the large transit depths in these lines, 
BJIO argued that oxygen and ionized carbon are preferentially 
heated to a temperature more than ten times higher than the 
temperature of hydrogen within a layer in the atmosphere lo- 
cated between ^2.25 R,, and the boundary of the Roche lobe 
at 2.9 Rp. 

Most analyses of the HD209458b UV absorption signatures 
to date have either been limited to first order deductions, such 
as the effective size of the absorbing obstacle (VM03,VM04) 
or been ba sed on com plicated first principle model s for the at- 
mosph ere (ITian et al. 2005 ; Garcia Munoz 2007 ; Murray-Clay| 
[eraL][2 009, Ben-Jaffel and Hosseini 2010). We beUeve fliat 
there is an important role for an intermediate class of models 
that satisfy some basic physical constraints, but parameterize 
aspects of the atmosphere that are difficult or impossible to 
predict accurately. For example, it is well established that the 
thermosphere of HD209458b should be composed primarily 
of H and H^, but the location of the transition from H2 to H is 
uncertain with d ifferent physical models making vastly differ- 
ent pr edictions ( Liang et al.|2003 Yelle|2004 Garcia Munoz 
er, the boundary between the atmosphere and 



2007) . Moreover, 



mterplanetary space is dependent upon the unknown strength 
of the stellar wind and the planetary and interplanetary mag- 
netic field and has yet to be modeled in a realistic fashion. 
Finally, although it seems well established that the tempera- 
ture is of the order 10,000 K, the precise value depends on the 
heating efficiencies, which have yet to be calculated, and ra- 
diative cooHng by minor species, which is not included in any 



of the models. Because of these uncertainties it is important 
to analyze the data in a way that makes clear what aspects of 
the atmosphere are constrained by the observations and which 
are not. 

In this paper we introduce a generic methodology that can 
be used to interpret UV transit light curves in stellar emission 
lines. We demonstrate this methodology by using a simple 
empirical model of the thermosphere to analyze the H I and 
O I ti-ansit depths of HD209458b summarized by BJIO. In 
Section |2] we introduce a model for calculating transit light 
curves for planets with extended atmospheres and discuss the 
basic features of the model thermosphere. In Sectionl3]we dis- 
cuss the H I transit depth measurements in detail andconfirm 
that they can be explained by absorption by atomic hydro- 
gen below the Roche lobe. Nevertheless, we also show that 
the core of the H Lya absorption line is optically thick up to 
the Roche lobe and that the atmosphere is evaporating. Fur- 
ther, we demonstrate that the disagreement between BJ08 and 
VM03 is due to differences in the treatment of the data and 
different definitions of the transit depth. 

In the rest of Section |3] we discuss the transits in the O I 
triplet lines and show that the empirical model thermosphere 
with a solar abundance of oxygen can be used to obtain transit 
depths that are statistically consistent with the observations. 
We also address the feasibility of the idea that energetic oxy- 
gen atoms are present in the thermosphere of HD209458b and 
present alternative ways to explain O I transit depths that are 
comparable to or larger than the full-width H I depth. We con- 
clude Section|3]by discussing a variety of different models for 
the thermosphere of HD209458b and use them to calculate 
transit depths. In Section fflwe discuss the feasibility of our 
assumptions and in Sectionpjwe summarize our findings and 
conclusions. 

2. METHOD 

2.1. Transit light curves 

In order to compare the observed UV transit depths with 
models of the thermosphere, we need to calculate the flux 
decrement of the stellar emission lines observed at Earth at 
different times t during the transit. The wavelength-integrated 
transit depth, Ts\ (f ), is given by: 



Tsxit)-- 



1 



djF, 



S(X-X')dX 



/(A,x)exp[-T(A,f,x)-T«M(A)]t/Arf (1) 



where ds is the distance to Earth, ^(A - A') is the instrument 
response function, 5\ corresponds to the spectral resolution 
of the instrument, /(A,x) is the specific intensity of the star 
expressed as a function of position x on the stellar disk of 
area A^/, r(A;f,x) is the optical depth due to the planet and its 
atmosphere along the line of sight (LOS) from the host star 
to the observer, t«m(A) is the optical depth of the interstellar 
medium (ISM) along the sightline to the star, and F,, is the 
out-of-transit flux within 5\ observed at Earth. We note that 
interstellar absorption does not affect the transit depth unless 
the ISM is optically thick. In that case the transmitted flux is 
so small that the observations cannot be used to calculate tran- 
sit depths. This affects the analysis of transit depths calculated 
from unresolved or partly resolved stellar emission lines. 

The optical depth T(A,f,x) is a complicated function of 
time and position, and in general equation ([T]i cannot be in- 
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tegrated analytically. In order to integrate the equation nu- 
merically, we adopted a Cartesian coordinate system with the 
origin at the center of the stellar disk. At time f, the coordi- 
nates of the center of the planet in this system are: 



d(t) = ^(t)^/l-cos^a(t)] 
b{t) = (^{t) cos (0 COS [a (f )] 



(2) 
(3) 



where d and b are the absolute values of the x and y coor- 
dinates of the planet, respectively, ^ is the distance between 
the centers of the star and the planet, / is the inclination, 
and a is the angle between the LOS to Earth and the line 
joining the centers of the planet and the star in the orbital 
plane. The y coordinate of the planet at the center of the tran- 
sit, be = ^ecos(0 expressed in terms of the stellar radius 
R^,, is known as the impact parameter. The angle a is re- 
lated to the true anomaly 6{t) of a given orbital position by 
a = \6c-9{t) \ where 6c is the true anomaly at the center of 
the transit. If the orbit is circular, obtaining a is straightfor- 
ward. For an eccentric orbit the distance ^ and true anomaly 
as a functio n of time can be calcu lated by solving Kepler's 
equation (e.g. |Koskinen et al.|2009[ ). 

2.2. Optical depth 

Assuming that the atmosphere is spherically symmetric, the 
optical depth t{v,z) at frequency v along a LOS with a tan- 
gent altitude z from the center of the planet is given by: 



t(v,z) = 2 



aJp,T)ns (r)rdr 



(4) 



where «, is the number density of the absorbing species s, Uy 
is the absorption cross section, p is pressure, T is temperature, 
and r is the distance from the center of the planet. We note that 
p and r depend on the altitude r and thus is a function of 
altitude. For individual absorption lines the cross section is 
given by: 



(p, ^) ■■ 



fo 



(j>^ (a,v) 



(5) 



where /„ is the oscillator strength and (f)^(a,v) is the Voigt 
function. The arguments of the Voigt function are given by: 



a = 



V = 



Al'D = 



2kT 



where is due to natural broadening, I'o is the laboratory fre- 
quency of the absorption line, is the mass of the absorbing 
species, and A^o is the Doppler broadening parameter. 

The relationship between the tangent altitude z and the x 
and y coordinates of a position on the stellar disk is given by: 



z\t,x,y)=[x-d(t)f+[y-b(t)f 



(6) 



This relationship can be used to calculate z for any point on 
the stellar disk and the LOS optical depth t(;/,z) can then be 
calculated by integrating equation (j4|. If the atmosphere is 
not spherically symmetric, the problem of obtaining the LOS 
optical depth as a function of position on the stellar disk is 
much more complicated. In that case, equation Q must be 
integrated separately for each latitude point, taking into ac- 
count the longitudinal variation of the atmosphere along the 
LOS. This technique may be relevant if, for instance, future 
observations of UV transit light curves indicate that the upper 
atmosphere deviates significantly from a spherical shape. 



As stated above, we integrated equation ([T]l numerically in 
Cartesian coordinates centered on the stellar disk. We used 
Simpson's rule to set up a sequence of one-dimensional inte- 
grals in the y direction that cover the disk of the star in the 
X direction. This technique is iterative and the algorithm in- 
creases the number of points on the stellar disk until a certain 
predetermined numerical precision is achieved [see Section 



4.6 of |Press et al.J ( |l992[ )]. The transit depths in the UV are 
of the order of 1-10 %, and thus we required numerical errors 
less than 0.1 % in all of our simulations. Given a model atmo- 
sphere, equation (|4]i was integrated by using the method de- 
scribed by Smith and Hunten ( 1990 ) (see Section 2.3 below). 
The optical depths were calculated for a fixed altitude grid 
of the model atmosphere, and the values were interpolated to 
the appropriate stellar disk coordinates by the iterative algo- 
rithm that was used to integrate equation ([T}. Voig t functions 
were evaluated numerically by using the method of |Humlicek| 
( il982| ). 

2.3. Model atmosphere 

Model atmospheres based on hydrodynamics and photo- 
chemistry are complicated and time-consuming to use. In 
addition, they are based on many uncertain assumptions that 
cannot always be changed easily to match with the observa- 
tions. In order to fit the data without the bias due to any indi- 
vidual model, we have constructed a simple empirical model 
of the upper atmosphere of HD209458b. This model is based 
on only a few free parameters that can be constrained by the 
generic features of the more complex models. Figure[T]shows 
the structure of the upper atmosphere and highlights some of 
the important transition altitudes of the model. The lower 
boundary of the thermosphere is located at the pressure of 
Pb = 0.1 /ibar, at an altitude z;, « 1.1 Rp. This choice is moti- 
vated by the fact that the stellar EUV radiation, which causes 
intense heating and ionization of the upper atmosphere, is al- 
most fully absorbed above the 1 /ibar level and a s a result the 
volume heating rate peaks near the 0. 1 //bar level (IKoskinen et 
al. 2010). At higher altitudes H2 is dissociated thermally and 
by dissociative photoionization. This is also true of the more 
complex molecules, and at z > 1.1 Rp th e atmospheric con- 
stituents appe ar as atoms or atomic ions ( |Yelle|2004{ |Garcia| 
|Munoz|2007l ). 

We assume that the lower atmosphere below Zb is opaque 
to the FUV radiation considered in this study, and that the 
mean effective temperature is T = 1300 K up to the 0.1 /ibar 
level. This temperature, together with the assumption of hy- 
drostatic equilibrium, is used to constrain the vertical extent 
of the lower atmosphere. We note that X rays and mid-UV ra- 
diation penetrate to the mbar levels, while the bulk of the vis- 
ible radiation penetrates down to the ^1 bar level or slightly 
above it. Thus these types of radiation do not have a sig- 
nificant effect on the thermosphere where the conditions are 
largely determined by EUV heating and ionization. 

Several different models suggest that the temperatu re of 
the thermosphere is T ^ 10,000 K (e.g., [ Ye lle 2004jjG^ 
cia Munoz||2007[ |Murray-Clay et al.|[2009[ [Koskinen et al. 



2010| ). Consequently, the pressure scale height is relatively 
large, reaching values of Hp > 7000 km. This implies that 
the thermal escape parameter A^. = r /Hp < 13, and thus sev- 
eral authors have suggested that the upper atmosphere of 
HD209458b is undergoing fast hydrodynamic escap^ (e.g. 

By fast hydrodynamic escape, we mean outflow that reaches supersonic 
velocity at some critical altitude. 
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Fig. 1. — An empirical model of the atmosphere of HD209458b. Several 
important transition altitudes are shown. Above z.b = 108 R,, (p = 10"' bar) 
molecules are dissociated, and atomic species of H, C, and O have been de- 
tected (VM03.VM04). The H+/H ionization front is located between Ze = 2- 
5 Rp and the ionization fronts of the other species may also be located at a 
similar altitude with H. The critical altitude where the outflow velocity be- 
comes supersonic in fast hydrodynamic escape is located at Zc ~ 3 R,,. The 
Hill radius is ~ 4.1 R,, and the distance of the Roche lobe from the cen- 
ter of the planet in the direction perpendicular to the LOS is Z/yx ~ 2.9 R,, 
{Gu et al . 2003 1. The vertical arrows show the different types of stellar radi- 
atron and their approximate penetration depths. In addition to hydrogen and 
helium the list of atmospheric constituents in the lower atmosphere includes 
the species that have been detected in the atmosphere of HD209458b (e.g. 
[Charbonneau et al..2002,.Knutson et al.,2008„Swain et al.,2009j . 



Lammer et al.||2 0031 iTTa neraLjlSoaSl [Garcia Munozl[2007l 
Murray-Clay et al. 2009 ). Despite this possibility we assume 
that the density profile is hydrostatic. This assumption is jus- 
tified because the density profile is approximately hydrostatic 
up to the critical level at Zc ~ 3 Rp even when the atmosphere 
is escaping at supersonic speed at higher altitudes. 

Based on the ratio of the vertical advection term to the 
pressure gradient term in the momentum equation, the out- 
flow velocity must be nearly equal to or faster than the speed 
of sound for the departure from hydrostatic equilibrium to 
be significant globally. Most models indicate that the mass 
loss rate from the atmosphere of HD209458b is limited to 
dM/dt = 10^-10*^ kg s"' by the available EUV energy. The 
pressure at Zc = 3 Rp in the empirical model shown in Fig- 
ure [T|is p PS 0.04 nbar, which implies that the vertical velocity 
at the critical level is = 1-10 km s"'. Due to the conser- 
vation of mass, the vertical velocity is inversely proportional 
to density and thus it decreases steeply with decreasing alti- 
tude. For example, at p = I nbar, which corresponds to the 
altitude of z = 1.7 Rp, the vertical velocity required to support 
mass loss is only v,- = 30-300 m s"'. Such a slow vertical ve- 
locity does not cause a significant de viation from hydrostatic 
equilibrium in the thermosphere (e.g., |Yelle|2004| l. 

Given a hydrostatic density profile, we integrated equa- 
tion (|4]) numerically. For the purposes of the numerical so- 
lution, we created a grid with an altitude spacing of 0.01 Hp 
for the thermosphere. The number density «, of species s at 
altitude corresponding to the grid point k is simply given 
by: 

(rk) = XsH (ri_i) exp[ At- (r^) - Xj {ru-i)] (7) 

where is the volume mixing ratio of species s, and n is the 
total number density. The thermal escape parameter, Aj-, is 



(8) 



given by: 

GMpink-i 
kTk-\rk 

where m^_i and 7]._i are the mean molecular weight and tem- 
perature, respectively, at a lower grid point k- 1. 

In order to match the observations we treated the number 
density at the lower boundary, «/,, the mean temperature of the 
thermosphere, T , and the mixing ratios as free parameters. 
In practice, we used pi, and T to constrain nt,, and assumed a 
multiple of the solar composition to fix x,. The resulting mix- 
ing ratios of H, He, C, O, and N were used to calculate the 
mean molecular weight m. The upper boundary of the model 
is at an altitude Ze above which the atmosphere is mostly ion- 
ized. The altitude Ze need not be the same for different species 
and generally it depends on the ionizing flux, photoabsorption 
cross sections, photochemical reaction rates, and the outflow 
velocity. We assumed that the the column densities of the 
neutral species are negligible above Ze and treated it as a free 
parameter for each species. 

We also assumed uniform mixing throughout the model at- 
mosphere and that diffusive separation does not take place 
in the thermosphere. This assumption can be justified be- 
cause the large pressure scale height in the atmosphere of 
HD209458b may lead to efficient turbulent mixing and im- 
plies that the exobase is located at a relatively high altitude. 
In general, diffusive separation occurs when the coefficient of 
turbulent (eddy) diffusion, K^, is smaller than the coefficients 
of molecular diffusion, D,, for different species. The coeffi- 
cient of eddy dif fusion can be estimated crudely as Kr ^ VrHp 
(e.g., Atreya 1986), where vv is the characteristic velocity of 
turbulence. The pressure scale height in the extended thermo- 
sphere of HD209458b is Hp > 10'' m and consequently it is 
possible that > 10^ m^ s"'. The coefficients o f molecular 
diffus ion reach similar values at < 20 nbar (Koskinen et al.| 
2010| l indicating that uniform mixing can take place up to the 
nbar levels. 

Drag forces between the escaping hydrogen and minor 
species can also lead to uniform mixing of the thermosphere. 
Whether or not a heavier atom moves along with hydrogen de- 
pends on its mass. If the mass of the atom is smaller than the 
crossover mass M^, the drag forces are strong enough to pre- 
vent diffusi ve separation. Appr oximately, the crossover mass 
is given by ( [Hunten et al.|1987] l: 

kTFn 



Mc = Mh + 



bgXH 



(9) 



where Mh is the mass of atomic hydrogen, Fh is the es- 
cape flux, and b is the binary diffusion coefficient for hy- 
drogen and the heavier species. Assuming a mass loss rate 
is dM/dt = 10^-10^ kg s~' and a binary diffusion coefficient 
of b ^ 1.6 X 10^^ m"' s"', we obtain the crossover mass of 
Me = 11-102. This implies that drag forces affect the mix- 
ing ratios of oxygen and carbon atoms, and could explain a 
significant abundance of these atoms in the thermosphere. 

3. RESULTS 

3.1. The distribution of hydrogen 

In this section we investigate the degree to which the ex- 
isting UV transit observations of HD209458b are consistent 



with the simple atmospheric model described in Section 2.3 
We begin by describing the fitting of the model to the tran 
sit depth measurements covering the full width of the stel 
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Fig. 2. — The model H Lya profile in the reference frame of the star The 
solid line is the unattenuated out-of-transit stellar flux and the dotted line 
shows the flux after the absoiption by the ISM has been applied. Two fea- 
tures are visible: a broad absorption band due to interstellar hydrogen (H I) 
between 1215.5-1215.9 A and a narrow absorption band due to interstellar 
deuterium (D I) centered at a laboratory wavelength of 1215.34 A. Both fea- 
tures are shifted according to the ISM parameters of |Wood et al.|j2005) . 

lar H Lya emission line profile in the wavelength range of 
[1212,1220] A. We note that VM04 obtained a Hne-integrated 
absorption depth of 5.3 ± 1.8 % in this wavelength range 
while BJIO obtained a higher depth of 6.6 ± 2.3 %. 

In order to calculate the depth of the transit, we used equa- 
tion (fTJ together with the absorption line parameters listed 
in Table [T] We simulated the stellar emissions by using the 
HLyg line profile of HD209458b calculated by [Wood et al] 
( 2005 p (hereafter W05). This line profile is based on a high 
resolution spectrum obtained with the STIS E140M grating 
and it was generated by reversing the absorption of the line 
profile by the ISM. We fitted a sum of three Gaussian distribu- 
tions to the pre-ISM line profile to represent the self -reversed 
core and the broad wings of the line. We then scaled the line- 
integrated flux to match with the observations of HD209458b 
analyzed by BJ07 after accounting for interstellar absorption. 

Most of the flux contained within the line profile is absorbed 
by the ISM, and thus it is important to make sure that the 
model fluxes are properly attenuated before the line-integrated 
transit depth is calculated. W05 measured the column density 
of hydrogen towards HD209458 and it is Nh = 2.3 x 10^^ m"^ 
They reported a Doppler broadening parameter along this line 
of sight of bo = 12 km s"'. This Doppler parameter is typ- 
ical of the local interstellar medium (LISM) where the tem- 
perature is thought to vary be tween Tjsm = 7,000-12,000 K 
(Landsman and Simon 1993 i. The column density, on the 
other hand, is relatively low and indicates that the sightline 
to HD209458 is in a low-density direction of the LISM. We 
calculated the extinction as a function of wavelength by using 
these parameters and adopting a Voigt function to simulate 
thermal and natural broadening in the ISM. Figure l2] illus- 
trates the effect of interstellar absorption on our model line 
profile. 

We obtained the observed full-width H I transit depth of 
6.6 % by adopting a mean temperature of T = 1 1,000 K and 
assuming that most of the hydrogen in the atmosphere of 
HD209458b is ionized above Ri 2.9 Rp. The rest of the 
free parameter values for this best-fitting empirical model, la- 
beled Ml, and the other models discussed in this paper are 



Fig. 3. — H Lya emission lines based on the Ml model (see Table[2j con- 
volved with an LSF corresponding to the spectral resolution of 0.08 A appro- 
priate for the STIS G140M grating. The solid line shows the out-of-transit 
flux and the dotted line shows the in-transit flux. The full-width transit depth 
based on these line profiles is 6.6 %. 



listed in Table |2] In order to illustrate the appearance of the 
transit in real observations, we convolved the data with the 
line-spread function (LSF) appropriate for the STIS G140M 
grating. Figure [3] shows the resulting synthetic stellar H Lya 
line in and out of transit. 

The location of the upper boundary in the Ml model agrees 
with models of the thermosphere, which suggest that the H/H"*" 
ionization front is located between = 1.5-5 ( Yelle 20041 
Garcia Munoz 2007). It is also close to the boundary 



of 



the Roche lobe in the direction perpendicular to the LOS at 
Zh± ^ 2.9 Rp. Further, the pressure at Z^, is p = 0.04 nbar, 
which is comparable to the pressure of the solar wind scaled 
to fl = 0.047 A U. Thus, if the planet lac ks a significant mag- 
netic field (e.g. |Griessmeie r et al. 2004 ), the upper boundary 
of the empirical model may be close to the ionopause. 

The optical depth profiles for the Ml model are shown in 
Figure |4| for the core of the line atA = 1215.67 A and for the 
wing of the line at A = 1215.2 A, which corresponds to Av « - 
1 16 km s"' in the stellar reference frame. We agree with BJIO 
that due to thermal and natural broadening, and large column 
densities along the lines of sight through the atmosphere, the 
optical depth in the wings of the line profile is significant even 
in the absence of actual bulk flows towards or away from the 
observer. The core of the line is optically thick up to and 
beyond the upper boundary while the wing is optically thick 
up to the altitude of z = 2 Rp. The latter altitude corresponds 
to a transit depth of ^5.4 %. 

Figure |5] shows transmission of the stellar H Lya emission 
during transit as a function of wavelength based on the Ml 
model. It also shows the transmission data points from Figure 
6 of BJ08. The Ml model is consistent with ^11 % absorp- 
tion within the core of the line profile. We calculated ~ 1-9 
for our fit to the data points. We note that BJIO o btained 
nearly the same value by using the DIVl model of [Garcia] 
|MunoZ| (2007) (hereafter GM07) after multiplying the density 
profiles oTthat model by a factor of 2/3. The results indicate 
that our fit to the data is as good as that obtained by BJIO. 
This is not surprising because in both the Ml empirical model 
and the DIVl model the mean temperature is T ss 10,000 K 
and the density profiles are approximately hydrostatic in the 



6 



T.T. Koskinen et al. 



TABLE 1 
Absorption line parameters" * 



Line 


Ao(A) 


/o 




Ao(s-') 


AAd(A) 


o-o(m2) 




HILyo 


1215.67 


4.164 X 


10-' 


6.265 X 10** 


0.087 


5.87 X 10- 


18 


DILyo 


1215.34 


4.165 X 


10-' 


6.270 X 10** 


0.062 


6.23 X 10- 


18 


OI 


1302.17 


4.887 X 


10-2 


3.204 X 10** 


0.023 


2.95 X 10- 


18 


01* 


1304.86 


4.877 X 


10-2 


1.911 X 10** 


0.024 


2.95 X 10- 


18 


01** 


1306.03 


4.873 X 


10-2 


6.352 X 10^ 


0.024 


2.95 X 10- 


18 



from |Morton1|T99T) and |Morton|(2003) 
Doppler EWHMs and central cross secti 



sections were calculated for T = 10,000 K 



TABLE 2 

Model parameters 



ID 


T(K) 


Z, (Rp) 


Pb ipbav) 


1% (m ^) 


m° 


Vr (km s ' ) 


Nm (m-2)'' 


H I (%) 


O I (%) 


(H 


Ml 


11,000 


2.92 


0.1 


6.6 X lO'* 


1 


0.0 


5.8 X 10" 


6.6 


4.3 


1.9 


M2 


11,000 


2.92 


0.1 


6.6 X 10'*" 


1 


10 


5.8 X 10^3 


6.6 


5.8 


1.9 


M3 


11,350 


2.92 


0.1 


6.4 X lO'* 


5 


0.0 


5.7 X 1023 


6.6 


5.0 


1.9 


M4 


12,000 


2.96 


0.1 


6.0 X 10'*" 


10 


0.0 


5.4 X 1023 


6.6 


5.5 


1.9 


M5 


15,950 


3.08 


0.1 


4.5 X 10"^ 


40 


0.0 


4.0 X 1023 


6.6 


7.1 


2.2 


M6 


11,000 


2.92 


0.001 


6.6 X lO"* 


1 


0.0 


6.1 X 102' 


2.3 


2.3 


8.9 


M7 


8250 


2.72 


1 


8.8 X 10" 


1 


0.0 


5.4 X 1024 


6.6 


3.9 


1.8 


KIO 


10,280'' 


2.92 


10 


6.6 X 10" 


1 


0.0 


6.8 X 1023 


5.2 


3.7 


3.1 


Y04 


12,170'' 


3.0 


200 


1.9 X 10^' 


1 


2.2'' 


9.9 X 1022 


2.9 


2.7 


7.3 



" Metallicity enhancement. In practice this factor is used to multiply the O/H and C/H ratios. 
^ The vertical column density of H I in the thermosphere. 

The chi-squared values are provided merely to aid comparison between different models. 
'' The model includes a self-consistent pressure-temperature profile. The quoted temperature is an average of the altitude levels. 
' The stated velocity applies to the upper boundary. 




1.5 2.0 
Altitude (Rp) 
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1216.0 



1216.5 



Fig. 4. — Optical depth of the H Lya line as a function of altitude for 
the best-fitting empirical model atmosphere of HD209458b. The solid line 
shows the optical depth at the central wavelength of Ao = 1215.67 A and the 
dotted line shows the optical depth at A = 1215.2 A. The latter wavelength 
con'esponds to Av = -116 km s-' in the stellar reference frame. 



Fig . 5 . — Transmission as a function of wavelength across the stellar H Lya 
emission line during the transit of HD209458b. The data points were obtained 
from Figure 6 of BJ08. The solid line shows our model fit to the data. The 
chi-squai'ed value for the fit is ~ 1.9. 



relevant altitude range of z = 1.08-3.0 Rp. 

Figure [6] shows the Ml model light curve as a function of 
time from the center of the transit compared with the observed 
data points from Figure 2(a) of BJ07. In order to calculate 
this light curve, we used the wavelength limits specified by 
BJ07 and BJ08. The transit depth at the center of the transit is 



^^6.6 %, which is slightly lower than that given by BJ07 but 
the difference is not statistically significant. The symmetry of 
the light curve and a good fit by our model indicate that there 
is no observational evidence that the extended thermosphere 
of HD209458b deviates significantly from a spherical shape. 
This is in contrast to the findings of VM03 who, based on 
the asymmetry of their light curve, argued that the planet is 
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n R h .... I .... I .... I .... -I 
-10000 -5000 5000 10000 

Time from transit (s) 

Fig. 6. — Model light curve showing the H I transit depth as a function 
of time from the center of the transit. The data points were obtained from 
Figure 2(a) of BJ07 and the solid line is the best-fitting model light curve. 
The appropriate wavelength limits specified by BJ07 and BJ08 were used to 
calculate the transit depths. 

followed by a curved cometary tail of escaping hydrogen. 

Next we explore whether the full width transit depth of 5- 
7 % is consistent with the larger depth of 8-15 % measured 
over a more limited wavelength range within the line profile 
(VM03,BJ07,BJ08). VM03 obtained the depth of 15 ± 4 % 
by comparing the ratio of the integrated flux within the wave- 
length intervals of [1215.15,1215.5] A and [1215.8,1216.1] 
A (the blue and red sides of the In region, respectively) 
to the flux within the intervals of [1214.4,1215.15] A and 
[1216.1,1216.8] A (the blue and red sides of the Out region, 
respectively) obtained during the transit. In their analysis, the 
interval [1215.5,1215.8] A was excluded because contamina- 
tion by the geocorona and sky background makes the flux 
measurements in this region unreliable. BJ07 redefined all 
of these wavelength intervals and obtained the transit depth 
of only 8.9 ± 2.1 %. In his analysis the 'contamination 
window' was extended to [1215.41,1215.84] A and the blue 
and red intervals were defined as [1214.83,1215.36] A and 
[1215.89,1216.43] A, respectively. 

Another significant difference between the analysis of 
VM03 and BJ07 is the treatment of the time variability of 
the H Lya emissions, either due to stellar variability or noise. 
BJ07 derived unperturbed H Lya profiles by merging the sub- 
spectra of a dense time series that was corrected for the transit 
trend. In the resulting line profiles (Figure 3 in BJ07), the 
time variability has been reduced to ~'l % level. VM03, on 
the other hand, attempted to account for time variability by 
calculating the transit depth based on the core-to-wing flux 
ratio instead of evaluating the transit depth based on the ratio 
of the flux obtained during transit to the flux obtained outside 
of the transit. This approach relies on the assumption that the 
transit depth in the wings of the line is not significant and that 
the variability in the wings is similar to the variability within 
the core of the line. 

Vidal-M adjar et al.| (j2008) (hereafter VM08) correctly 
pointed out that a comparison between the two approaches 
is not fair because BJ07 used different wavelength limits and 
instead of evaluating the change in the core-to-wing (In/Out) 
flux ratio, he calculated the true transit depth i.e., the flux 



decrement based on equation (iTli. However, we believe that 
the transit depths calculated by vM03 are not necessarily re- 
lated to the size of the absorbing cloud while those calculated 
by BJ07 and BJ08 are. By using the Ml model and the VM03 
and BJ07 wavelength limits, we obtain absorption depths of 
2.4 % and 1.3 %, respectively, based on the In/Out flux ratio 
and true transit depths of 7.5 % and 6.6 %, respectively. Our 
model is a good fit to the transmission curve of BJ08 and the 
behavior of the In/Out ratio is roughly in line with the line 
profiles shown in Figure 3(a) of BJ07 and yet the transit depth 
based on the In/Out flux ratio is small. We note, however, that 
the Ml model overflows the Roche lobe and thus we agree 
with VM08 that the atmosphere is evaporating. 

3.2. The distribution of oxygen 

The large transit depths in the O I and C II lines reported 
by VM04 and BJIO can potentially be used to constrain the 
degree of eddy mixing and the composition of the upper at- 
mosphere. The O I transit depth, if accurately measured, is 
a good indicator of the proportion and mixing of heavy el- 
ements in the thermosphere. The interpretation of the C II 
transit depth, on the other hand, is much more complicated 
because of the uncertainties related to the ionization profiles, 
photochemistry and, if C II is escaping the thermosphere, the 
interaction of the atmosphere with the stellar wind or mag- 
netospheric plasmas. Statistically, the observed C II transit 
depth is also the least significant detection out of the three 
currently available UV transit depth measurements. For these 
reasons, we limit our discussion and analysis to the interpre- 
tation of the O I transit depth. 

VM04 reported an integrated transit depfli of 12.8 ± 4.5 % 
in the wavelength range of [1300,1310] A, which covers the 
O I triplet emission lines listed in Table [T] BJIO reported a 
revised and slightly shallower transit depth of 10.5 ± 4.4 %. 
These transit depths are only 2.4 and 1.9 a, respectively, away 
from the FUV continuum transit depth of ^2 %. The data 
were obtained by using the STIS G140L grating with a spec- 
tral resolution of ~'2.5 A, which does not allow for the in- 
dividual lines of the triplet to be resolved. The lines were 
resolved, however, in an earlier spectrum obtained by VM03 
with the E140M echelle grating. Although this spectrum was 
not used to analyze the transit depths, it is a good reference 
point for the expected line profiles and can be used to con- 
strain the properties of the ISM. 

We modeled the emission spectrum of the O I triplet by us- 
ing the full disk solar line profile parameters derived by Glad^ 
stone ( 1992) and fitted the line-integrated fluxes to the obser- 
vations. As a result we obtained surface fluxes of 3.0, 3.2, and 
3.4 W m-2 for the O I, O I*, and OP* fines of HD209458, re- 
spectively. In order to correct the model fluxes for absorption 
by interstellar oxygen, we assumed that the O I/H I ratio in the 
LISM is -3 X 10"^ (Moos et al. 2002) and that afl of the in- 
terstellar oxygen is in its ground state. The latter assumption 
is consistent with the ratios of the individual emission lines in 
the high-resolution E140M spectrum. 

BJIO used the E140M data points to deduce fluxes for the 
O I, O I*, and O I** lines that scale as 1:1.5:1.17. The so- 
lar line-integrated fluxes of these lines scale as 1:1.06:1.16. 
After applying interstellar absorption to the solar line pro- 
files, we obtained model fluxes that scale as 1:1.43:1.55. The 
close agreement between the model and observed O I/O I* 
flux ratios confirms that most of the interstellar oxygen is in 
its ground state and that our parameters for the LISM are re- 
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liable. The observed O I**/0 I flux ratio, on the other hand, 
is slightly lower than the model flux ratio. This cannot be due 
to absorption in the ISM because a population inversion in 
the excited levels is unlikely. Thus the disagreement, if real, 
arises from the intrinsic properties of HD209458 and implies 
that the stellar O I triplet flux ratios differ slightly from the 
corresponding solar flux ratios. We note, however, that the 
E140M spectrum is noisy and it is not clear if the difference 
is real or if it arises from uncertainties in the data. 

The flux in the core of the O I ground state line is ab- 
sorbed by the ISM, and a significant O I transit depth is pos- 
sible only if the populations of the excited states of oxygen in 
the thermosphere of HD209458b are sufficient. We estimated 
the relative population of the ground state and the excited 
states by assuming that the thermosphere is in LT E and by 
using the electron ic partition function calculated by |Coloniia| 
[and Capitel lT(^9). At a temperature of T = 10,000 K we 
obtained population ratios of 5.2:3.0:1 for the ground state, 
lower excited level, and upper excited level, respectively. 
These ratios agree closely with the statistical weights of the 
levels that give rise to the triplet. 

The Ml model presented in Table |2] yields an O I transit 
depth of 4.3 %. This is similar to the transit depth calcu- 
lated by BJIO who used the scaled density profiles from the 
DIVl model of GM07. Again, the agreement between these 
models is not surprising because the temperature and density 
profiles are similar and both models assume solar O/H ratios. 
We note that the O I transit depth based on the Ml model is 
within ^1.4(7 from the observed depth. Assuming Gaussian 
statistics, the probability that the true transit depth is between 
1.7-6.1 % (2-1ct) is ~14 %. Due to the relatively large un- 
certainty in the data points, this probability is significant. In 
other words, the deviation of the model transit depth from the 
observed value is not statistically significant. Better observa- 
tions with higher S/N and spectral resolution are required to 
constrain the transit depth further Our analysis shows that the 
STIS G140M grating, rather than the G140L grating, could 
already be used to obtain better data. 

BJIO argued that the O I transit depth of 10.5 % is possible 
if oxygen is preferentially heated to Tqi ^ 100,000 K (8.6 eV) 
within a layer located between z = 2.25 and Zh± = 2.9 Rp. 
It is true that hot oxygen populations have been observed 



1993^ Shematovich 



in the thermospheres of terrestr i al planets (e .g. Cotton et al. 

eraL][2006| |Nagy et al.p981 |l but thes'e 
led mostly by dissociative recombination 



populations are formed mostly by 
of O2, which should not be important on HD209458b. Other 
possible mechanisms for heating oxygen atoms include colli- 
sional processes that, as we demonstrate in the next section, 
would also imply that oxygen is mostly ionized. Thus it is 
not clear if a neutral cloud at these temperatures is physically 
possible. 

As we already pointed out, the O I transit depth of 4.3 % 
based on the Ml model is consistent with the observations. 
However, the measurements also allow for a larger transit 
depth and it is thus worth exploring if transit depths close to 
10.5 % can be explained by changing the parameters of the 
model atmosphere, the host star, or the interstellar medium. 
In the following we investigate several different possibilities 
that can give rise to larger O I transit depths without the need 
to introduce energetic oxygen atoms in the thermosphere. 

Theoretically, the O I transit depth is sensitive to interstel- 
lar absorption. In order to test if changes in the abundance of 
interstellar oxygen can increase the calculated transit depth, 
we assumed that the column density of interstellar oxygen to- 



wards HD209458b is negligible. As a result, the calculated 
transit depth increased only by fractions of a percent. We con- 
clude that, contrary to H I, the full-width O I transit depth of 
HD209458b is not very sensitive to interstellar absorption. In 
addition, the E140M spectra indicate that our ISM parame- 
ters are accurate and the assumption of a negligible oxygen 
column would not be reasonable anyway. Consequently, we 
must explore the properties of the host star and the planet in 
order to explain the possibility of large O I transit depths. 

Transit depths in the visible and the IR wavelengths are sen- 
sitive to stellar activity and can appear larger if a significant 
portion of the stellar disk is covered by starspots - even if 
these spots do not actually fall in the path of the transit (e.g.. 



Czesla et al. 2009 1. This is because the effective emitting area 
of the star appears smaller In order to crudely estimate the 
magnitude of this effect, we defined an effective fractional 
area, Ag, of the quiet regions as the ratio of the area that would 
be required to emit zero flux to account for the diminished lu- 
minosity during quiet periods to the area of the stellar disk. In 
terms of A^^, the observed transit depth, dr, can be written as: 



i-A„ 



(10) 



where Rp is the radius of the transiting planet, and is the 
radius of the star. We note that Aq would have to be ^60 % 
for the O I transit depth to increase from 4.3 % to 10.5 %. 
Because the definition of Ag implies that the actual fractional 
area of the quiet regions would have to be much larger than 
60 %, this is clearly not realistic. In fact, the O I tripl et flux 
of the Sun varies by ^30 % during the solar cycle (Glad- 
p92| 



stone 



most by 



and thus we would expect an enhancement at 
percentage points (even this is probably overes- 



timated because the relationship between the flux variability 
and is complicated). This is not statistically significant and 
we conclude that the interpretation of the low-resolution tran- 
sit depths is not affected by the distribution of active or quiet 
regions. 

The transit depth measurements are also affected by stel- 
lar limb darkening or limb brightening. According to the 
solar disk ob servations of the O I triplet fluxes reported by 
|Roussel-Dupre^ (1985) , the fluxes are roughly constant be- 
tween /j,(= cos^) ~ 1--0.2, with only negligible limb darken- 
ing. At /i < 0.2 the fluxes are brighter than within the center 
of the disk by a factor of ^1.1. If the atmosphere of the tran- 
siting planet does not cover the bright ring surrounding the 
stellar disk, limb brightening actually makes the transit depth 
appear shallower We simulated limb brightening by dividing 
the stellar disk into a central disk region and a bright ring, 
and neglected limb darkening within the central region. As 
a result, we found that limb brightening does not lead to a 
detectable change in the low-resolution O I transit depth. 

Several authors have attempted to explain the observed O I 
transit depth by referring to the velocity dispersion arising 
from radial outflows. VM04 invoked this explanation and ar- 
gued that hydrodynamically escaping oxygen atoms moving 
at sonic or supersonic velocities could give rise to a transit 
depth of 13 %. Later, Tian et al. (T005) argued that the veloc- 
ity dispersion arising from similar supersonic outflows could 
also explain the large H I transit depth reported by VM03. 
They assumed that the Doppler broadening parameter is given 
by: 

l2kT 



(11) 
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where rj is the velocity dispersion. Assuming that 
77 = 10 km s"', and using the Ml model as a base, we obtain 
H I and O I transit depths of 6.7 and 7.7 %, respectively. Chi- 
squared for the H I transmission fit in this case is ^ 2.1, 
which indicates that the H I transit depth is not significantly 
affected by additional velocity dispersion. This is because a 
large portion of the H Lya line is akeady saturated. On sur- 
face, then, velocity dispersion appears to be a promising way 
to explain the large O I and C II transit depths. 
Equation 
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however, is appropriate for simulating 
broadening due to microturbulence. It is only accurate if the 
scale of the turbulence is much smaller than the mean free 
path. Further, what really matters is not radial velocity but 
the projected velocity along the LOS. Because most of the 
absorption during transit is due to the layers near the termina- 
tor, radial velocity is not a good measure of the true velocity 
dispersion. In order to estimate the velocity dispersion prop- 
erly, we calculated the appropriate velocities for each layer 
along each LOS and used these velocities in integrating equa- 
tion Q. For each layer k along a LOS, the parallel velocity 
v^ll is given by; 



Vk\\{z) = Vr{,n)co5[- 



-cos 



rk 



(12) 



where is the radial velocity, and z is the constant tangent 
altitude of the LOS. Assuming a uniform radial velocity of 
\'r ^ 10 km s"' for the Ml model, the maximum LOS veloc- 
ity, obtained at the upper boundary of the model along a sight- 
line through the lower boundary is Vao\\{Zb) ^ ±9.2 km s~'. 
For a LOS at z = 2 Rp, the corresponding velocity is only 

The M2 model (see Table [2| includes a uniform radial ve- 
locity of Vr = 10 km s"' . Figureplshows the transmission pro- 
file of the O I* line for this model together with a correspond- 
ing profile for the Ml model and another profile calculated 
by using equation ( [TT| . The absorption profiles differ signif- 
icantly. By using equation ( [TT| i we obtain a much broader 
absorption profile with a sligEtly smaller central cross sec- 
tion compared to the Ml model. If the velocity dispersion is 
treated properly, as it is in the M2 model, the absorption pro- 
file exhibits two symmetric peaks around the central wave- 
length of the O I* line at = 1304.86 A. The resulting H I 
and O I transit depths based on the M2 model are 6.6 and 
5.8 %, respectively. 

A radial velocity of = 10 km s"' is not enough to explain 
an O I transit depth of 10.5 % but it can explain O I tran- 
sit depths that are comparable to the H I transit depth even 
when the abundance of oxygen is significantly lower than 
that of hydrogen. It is not realistic, however, to assume that 
V,. = 10 km s"' in the lower thermosphere. Hydrodynamic 
models of the thermosphere indicate that the base of the plan- 
etary wind is located near p = \ nbar, which corresponds to 
an altitude of z = 1.7 in the Ml model. Assuming that the 
outflow is effective above this altitude, we obtain an O I tran- 
sit depth of only 5.5 %. Consequently, velocity dispersion can 
only affect the transit depth significantly if the optical depth 
above the critical level of Zc = 3 Rp is significant. 

It should be noted that radial velocities higher than 
Vr = 10 km s"' do not necessarily lead to larger transit depth 
because the absorption is limited by the fact that the central 
cross section in equation (|5]l is proportional to /S.vd- In fact, 
if we introduce a radial velocity of Vr = 60 km s"' to the Ml 
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Fig. 7. — Transmission as a function of wavelengtli in the O I* line. Tlie 
solid and dotted lines show transmission based on the Ml and M2 models, 
respectively (see Table|2j. The M2 model includes a uniform radial velocity 
of V,- = 10 km s"' and thus the absolution line profile exhibits two symmet- 
ric peaks around the central wavelength of Ao = 1304.86 A. The dashed line 
shows transmission for a model that includes the same radial velocity disper- 
sion but by using equation for microturbulence. 



model, we obtain an O I transit depth of only 5.2 %. Within 
the velocity range of vv = 10-60 km s"' , we obtain a maximum 
O I transit depth of 7.2 % by assuming that vv = 30 km s"'. 
Thus fast hydrodynamic escape of neutral oxygen from the 
thermosphere alone cannot explain an O I transit depth of 
10.5 %. 

The next logical step is to investigate if the O/H and C/H 
ratios in the thermosphere of the planet are actually higher 
than the solar values. BJIO did not consider this possibility 
because the density profiles of their best-fitting model are al- 
ready enhanced with respect to the Y04 model. With respect 
to the DIVl model of GM07, however, their density profiles 
are actually depleted and further, the O/H and C/H ratios of 
their model are not enhanced but solar. We generated new em- 
pirical models M3, M4, and M5 (see Table [2]) by increasing 
the O/H and C/H ratios of the Ml model by factors of 5, 10, 
and 40, respectively. Initially, this caused both the H I and O I 
transit depths to decrease because the higher mean molecu- 
lar weight of the thermosphere led to a shorter pressure scale 
height. In fact, we found that it is not possible to increase 
the abundance of the heavier atoms without also increasing 
the mean temperature of thermosphere to recover the fit to the 
H I transit depth. In order to do so, we increased the mean 
temperatures of the M4 and M5 models to T w 12,000 K and 
16,000 K, respectively. The M3 model did not require a sig- 
nificantly increased temperature. 

The mean temperature of the thermosphere can be con- 
strained by two additional means. First, the temperature of the 
neutral atoms is limited by thermal ionization. In LTE condi- 
tions the thermal ionization balance is determined by the Saha 
equation, which implies that atomic hydrogen is fully ionized 
at temperatures higher than T = 7000-8000 K. In non-LTE 
conditions of the thermosphere the ion/neutral fractions re- 
sulting from thermal ionization are determined by a balance 
between impact ionization by thermalized electrons and radia- 
tive recombination. The reaction rate coefficie nts of these pro- 
cesses, respectively, for atomic hydrogen are ( |Voronov|1997[ 
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Storey and Huminer|1995 1: 



7/0.39 

k = 2.91 X 10""^ exp(-U) s'^ 

U + 0.232 ^ 



aA = 4.2 X lO-'^lOV?;)' 



(13) 



(14) 



where is the electron temperature and U = Ip/Eg is the ra- 
tio of the ionization potential Ip to the energy of the electrons 
Ee « kTg. The balance of these rate coefficients implies that 
[H+] / [H] > 10 for T, > 25,000 K. Based on the impact 
ionization ra te of|Voronov|(| 19971 and the radiative recombi- 
nation rate of Schunk and Nagy (2000 1, the ionization fraction 
of atomic oxygen i.e., [O^J / [OJ > 100 in the same tempera- 
ture range. 

Second, as BJIO also points out, higher temperatures lead 
to a worsening fit to the H 1 transmission data in Figure |5] As 
the temperature approaches T = 15,000 K, > 2 for the fit. 
This is not a very strong constraint at all, of course, but it indi- 
cates a general trend that disfavors temperatures much higher 
than T = 15,000 K. Indirectly, the need to retain the fit to the 
H 1 transmission data and these temperature limits imply that 
models with O/H and C/H ratios much higher than 40 times 
solar ratios are not realistic, and that enhanced oxygen abun- 
dance can increase the O 1 transit depth at most to 7.1 %. We 
also experimented with a combination of supersolar oxygen 
abundance and hydrodynamic escape by introducing a uni- 
form radial velocity of = 30 km s"' to the M5 model. As 
a result, we obtained an O I transit depth of 10.8 %, which is 
actually close to 10.5 %. Despite being a good fit to the obser- 
vations, there are several good reasons for why this model is 
unrealistic. The O/H ratio of the M5 model is 40 times higher 
than the solar ratio and v,- = 30 km s~' below the critical level. 
Matching the observations is not enough prove a theory - the 
theory also has to make physical sense. 

It is possible that, contrary to our assumptions, there are 
optically thick clouds of neutral hydrogen and oxygen above 
the upper boundary of the Ml model moving at large LOS ve- 
locities. These clouds can be formed, for instance, by neutral 
atoms that escape the atmosphere and react with the stellar 



wind ( jHolstrom et al.l2008t |Ekenback e t al.|20101 l. Although 
the density profile of the empirical model is misleading at 
high altitudes, we extended the Ml model thermosphere to 
z = 10 Rp. As a result, we obtained H 1 and O 1 transit depths 
of 13.1 and 12.2 %, respectively. In this case, = 13.4 for 
the H 1 transmission fit. In a slightly more realistic case, we 
assumed that only 10 % of the atoms are neutral between 
and z= 10 Rp. This led to H 1 and 1 transit depths of 7.8 and 
6.3 %, respectively, and = 2.1. We note that external neu- 
tral atoms can affect the transit depth but we also agree with 
BJIO that the H 1 transmission data points are not consistent 
with optically thick populations of hydrogen moving at large 
velocities. Detailed models of the chemical and dynamical 
interaction of the escaping atmosphere with the stellar wind 
outside the Roche lobe are required to study the distribution 
of external neutral atoms. 

3.3. Models for the thermosphere of HD209458b 

With the exception of GM07, most of the modeling efforts 
so far have concentrated on simulating the distribution of hy- 
drogen (and sometimes helium) in the escaping atmosphere 
of HD209458b. These models are justified by the assumption 
that the overall characteristics of the thermosphere are not af- 
fected by small abundances of heavier neutral atoms or ions. 



In addition to the empirical model presented in Section |2.3[ 
we have used two models with different characteristics and 
limiting assumptions to calculate transit depths. These models 
are the one-dimensional outflow model of Yelle (2004) (Y04) 
and a t hermospheric circulation model developed by KoskH 
nen et al. (2010) (KIO). The photochemistry in these models 



is based on a composition of hydrogen and helium only. In 
order to roughly estimate the O I transit depth, we added a 
solar proportion of oxygen to these models and assumed an 
oxygen ionization fraction equal to that of hydrogen. 

The KlOa model is different from the existing outflow 
models in that it is three-dimensional and can account for 
horizontal flows that are important in the thermosphere of 
HD209458b. On the other hand, the model assumes that the 
thermosphere is in hydrostatic equilibrium and thus it can- 
not properly account for vertical flows that lead to a deviation 
from hydrostatic conditions. Also, the photochemical scheme 
of the model is not valid in regions where H"^ is the dominant 
species. To account for these limitations, we placed the up- 
per boundary of the transit depth calculations to z = 2.9 Rp 
in line with the arguments we have made for the validity of 
the empirical model. In order to calculate the optical depth as 
a function of altitude, we averaged the density and tempera- 
ture profiles of the model globally to obtain one-dimensional 
vertical profiles. The resulting H I and O I transit depths are 
5.2 and 3.7 %, respectively, and for the H I transmission fit 
we obtain ^ 3.1 (see Table pi). These values are in rough 
agreement with our earlier calculations and the observations. 

The one-dimensional outflow model developed by Y04 
solves the vertical component of the Navier-Stokes equations 
and includes a self-consistent treatment of much of the im- 
portant physics in the thermosphere of HD209458b. The alti- 
tude range of the model is limited to z < 3 Rp because above 
this altitude stellar gravity, radiation pressure, and interac- 
tion of the escaping atmosphere with the stellar wind affect 
the solution in ways that cannot be accounted for by sim- 
ple one-dimensional, hydrodynamic models. In order to en- 
sure consistency with kinetic theory, the velocity at the up- 
per boundary of the model was required to be equal to the 
weighted mean of the velocity distribution functions of differ- 
ent species. This is not the same as setting the velocity to be 
equal to Jeans effusion velocity because drifting Maxwellians 
were used to calculate the boundary condition in an iterative 
fashion. 

The density and temperature profiles of the Y04 model 
yield H I and O I transit depths of 2.9 and 2.7 %, respectively 
(see Table |2|. While the full-width H I transit depth is within 
1 .6(7 of the observed value, x^ for the HI transmission curve 
is larger than the one obtained for the empirical model or the 
KlOa model. The reason for the slightly poorer fit lies in the 
details of the density profiles. The H2/H dissociation front of 
the Y04 model is located at the pressure of p = 0.5 nbar at 
the altitude of z = 1.1 Rp. This means that the model is less 
extended than the empirical model below the 0.5 nbar level 
because the mean molecular weight of the gas is larger and 
the thermosphere is cooled effectively by infrared emissions 
from H3 ions. Also, the ionization front of the model is lo- 
cated at a relatively low altitude of « 1.7 Rp. 

The extent and temperature of the thermosphere depend 
strongly on the location of the H2/H dissociation front that 
effectively determines the number density and temperature at 
the lower boundary of the empirical model. In order to illus- 
trate this point, we changed the pressure at the lower bound- 
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ary of the Ml model from pi, = Q.l /ibar to = 1 nbar (see 
model M6 in Table |2]i. As a result we obtained reduced H I 
and O I trans it depths of 2.3 and 2.2 % , respectively. This 
explains why |Murray-Clay et al.| ( |2009] l obtained a relatively 
low H 1 transit depth of ^3 % by usmg a one-dimensional out- 
flow model with a composition based on H and H"^ only. The 
pressure and temperature at the lower boundary of the model 
are only p « 30 nbar and 7], = 1000 K, respectively. Thus the 
model ignores the critical region of the atmosphere that con- 
tributes significantly to the optical depth. If the pressure of 
their lower boundary was higher, they could have obtained a 
larger transit depth. 

In the KlOa model the H2/H dissociation front is located 
at the pressure of p = 0.1 /ibar due to a combination of ther- 
mal dissociation and horizontal mixing. In the DlVl model 
of GM07 it is located at p = 1.5 /ibar because, contrary to 
the other models, it includes ion and neutral photochemistry 
with carbon, oxygen, and nitrogen species. Thus H2 is also 
dissociated in reactions wit h atomic oxygen an d the OH radi- 
cal. We note, however, that Liang et al. ( 2003 1 also modeled 
the neutral photochemistry and found that H2 was not entirely 
converted into H at these pressure levels. Motivated by the 
DlVl model, we moved the lower boundary of the Ml model 
to pi,= I /ibar and obtained H 1 and O 1 transit depths of 6.6 
and 3.9 %, respectively (see model M7 in Table |2]). In this 
case, we adjusted the mean temperature of the thermosphere 
to r = 8250 K in order to retain the fit to the H 1 transit depth. 
In terms of optical depth, the M7 model is almost indistin- 
guishable from the Ml model. 

Our exploration together with the results of B J 1 reveal that 
models of the thermosphere can be used to match the observed 
H 1 absorption depths as long as the location of the H2/H dis- 
sociation front, the temperature and the location of the upper 
boundary agree with the values derived from our empirical 
model. These findings provide strong justification for our ap- 
proach of using a simple empirical model to fit the data. We 
would not have been able to identify these important parame- 
ters that affect the optical depth profiles in the thermosphere if 
we had allowed one of the more complicated models to cloud 
our judgement. Once known, these parameters can help to 
guide the hydrodynamic and photochemical models towards 
solutions that agree with the observed transit depths. 

4. DISCUSSION 

In this section we assess the feasibility of our basic assump- 
tions and results. We find that it is easy to generate models 
that agree with the observed H I transit depths. According to 
the best-fitting empirical models Ml and M7 (see Table |2]i, 
the H2/H dissociation front is located between p = 0.1-1 /ibar 
(z w 1.1 Rp) and the mean temperature of the thermosphere is 
between T = 8000-1 1,000 K. We note that due to thermal ion- 
ization, the temperature of the neutral thermosphere cannot be 
much higher than T « 15,000 K. The Ml model atmosphere 
is optically thick up to z « 2 Rp in the wings of the H Lya line 
and at least up to the Roche lobe in the core of the line. Above 
our upper boundary at = 2.9 Ry„ however, the atmosphere 
is mostly ionized. 

Our results are sensitive to the pressure level of the H2/H 
dissociation front. Additional observations are required to 
constrain the vertical column density of H I in the t hermo- 
sphere. Suc h observations were recently obtained by France 
et al. ( 2010| l who attempted to observe the auroral and day- 
glow emissions of H2 from HD209458b. They did not detect 
the emissions and argued that this impHes that the vertical col- 



umn density of H I must be at least Nhi = 3 x 10"^ 
prevent stellar FUV radiation from exciting H2 in the lower 
thermosphere. Table ^includes substellar column densities of 
H I for the models discussed in this paper. According to this 
additional constraint, it is possible that the H2/H dissociation 
front is deeper than the 0. 1 /ibar level. In that case, M7 is the 
best-fitting model. 

The detection or non-detection of infrared emissions from 
H3 can also be used to const r ain th e properties of the ther- 
mosphere. [Shkolnik et al.| ( |2006 ) attempted to observe 
these emissions from a sample that included the transiting 
planet GJ436b and several non-transiting EGPs by using the 
CSHELL spectrograph at the NASA Infrared Telescope Facil- 
ity (NIRTF). They reported upper limits for the fluxes that are 
often several orders of magnitude highe r than the predicted 
fluxes for close-in EGPs. More recently, |Swain et al. ( |2010) l 
obtained a near-IR secondary eclipse spectrum of HD 1 89733b 
by using the SpeX instrument at NIRTF. The lack of a 3.9 /im 
feature in this spectrum may imply that H3 emissions are not 
important on this particular planet. We note that a lack of 
emissions implies that there is not enough H2 to form H3 
or that H3 is destroyed in reactions with oxygen or carbon 
species. 

Modeling the chemistry and dynamics of the lower ther- 
mosphere is a complicated problem. For instance, GM07 has 
shown that the complex photochemistry that includes oxygen 
and carbon species affects the location of the H2/H dissocia- 
tion front and thus the extent and temperature of the thermo- 
sphere. H2 can also be dissociated by high-energy electrons 
that penetrate below the 1 /ibar level and have not been in- 
cluded in the existing photochemical models. The composi- 
tion of the atmosphere is also affected by dynamics. Because 
the upper atmosphere is strongly ionized, the conductivities 
are likely to be high. As a result, the dynamics of the ther- 
mosphere is affected by strong electric currents that have not 
been included in the existing models either 

The results are not as sensitive to the location of the upper 
boundary as they are to the location of the lower boundary. 
We note, however, that the location of the upper boundary at 
Ze = 2.7-2.9 Rp in both our work and that of BJIO is remark- 
ably consistent. It is tempting to suggest that this is not a 
coincidence. After all, the upper boundary is close to both the 
boundary of the Roche lobe and the ionopause of a weakly 
magnetized planet. It is possible that, for some reason, the 
density of the neutral atoms decreases sharply above this alti- 
tude. We believe that multidimensional plasma models, either 
kinetic or MHD, are needed to study the density profiles of 
different species above the upp er boundary. 

Based on an earlier study b ylHolstrom et al. (2008 1, |Eken-| 
[back et aL] ( |2010| ) (hereafter ElO) studied the formation of 
energetic neutral atoms (ENAs) outside the Roche lobe of 
HD209458b. They pointed out that the ENAs are created by 
charge exchange between hydrogen atoms escaping the atmo- 
sphere and the protons of the stellar wind. They argued that 
absorption by ENAs can explain the H I transmission curve re- 
ported by VM03. In principle, this is a good idea but some of 
the details are unclear We are not convinced that the bound- 
ary conditions of the ElO model are consistent with realistic 
models of the thermosphere. Also, the ElO model was not 
used to fit the BJ08 transmission data. 

The inner boundary of the ElO model is located at z = 2.8 Rp 
and the optical depth of the atmosphere in the H Lya line 
is assumed to be small inside the inner boundary sphere. 
At the same time the number density at the inner bound- 
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ary is « = 4 x lO'-' m ^ because a mass loss rate of 
dM/dt ^ 1.6 X 10" g s"' is required to ensure a sufficient 
production rate of the ENAs. We note that the number den- 
sity of the Ml model at z = 2.9 Rp is n = 2.6 x 10'^ m"^ and 
that the density profile of this model already leads to a full- 
width H 1 transit depth of 6.6 %. Thus absorption by ENAs is 
not required to fit the data. On the contrary, too much addi- 
tional absorption may actually lead to a deteriorating fit to the 
detailed H 1 transmission data. 

The mass loss rate of dM/dt ^ 1.6 x 10" g s"' is in line 
with the GM07 model. BJIO showed that this model can 
be used to fully explain the observed absorption. They also 
pointed out that an optically thick external population of neu- 
tral hydrogen is possible but the LOS velocities must be slow. 
In this case, additional absorption is constrained to the core of 
the line, which is fully absorbed by the ISM and the transmis- 
sion in the wings is not significantly affected. Our analysis 
supports this conclusion. New models that treat the thermo- 
sphere and the external region self-consistently are required 
to study the true distribution of external neutral atoms. 

The 1 transit depth is less precise than the H 1 transit depth 
but we are able to explain the observations with our empirical 
model and 0/H ratios ranging from solar to 40 times solar, 
although the latter value implies an uncomfortably high tem- 
perature. The large O 1 transit depth indicates that the ther- 
mosphere is uniformly mixed up to the Roche lobe although 
there is a marginal possibility that this is not the case because 
the FUV continuum transit depth of 2 % is within the 1 .96a 
confidence interval of the detected depth. If diffusive separa- 
tion takes place, the O I transit depth cannot be much higher 
than the continuum transit depth. 

Based on the Ml model, we obtain a full-width O 1 transit 
depth of 4.3 %, which is statistically consistent with the ob- 
served absorption. In this case, oxygen overflows the Roche 
lobe with hydrogen. Assuming radial velocities of vv =1- 
10 km s~' at = 2.9 Ry„ the total mass loss rate from the 
Ml model is dM/dt ^ 4^0 x 10^ kg s"' and this leads to a 
total escape flux of Fq = 1.8-18 x 10"^' s~' for oxygen. We 
have also shown that the velocity dispersion due to radial out- 
flow, supersolar oxygen abundance and external populations 
of neutral oxygen can explain O 1 transit depths that are com- 
parable to the full-width H 1 depth. We have not, however, 
been able to feasible mechanisms that give rise to 1 

transit depths that are significantly higher than the H 1 depth. 

We believe that more precise measurements are required to 
further constrain the 1 transit depth. We have generated syn- 
thetic data by simulating the response function of the STIS 
G140M grating to show that, by assuming the Ml model to 
be accurate, we would observe a transit depth of 4.3 ± 1 % 
for HD209458b. Thus for transit depths higher than -5 %, a 
confident detection in the upper atmosphere can easily be ob- 
tained. In addition, the spectral resolution of the G140M grat- 
ing is high enough to resolve the individual emission lines of 
the O 1 triplet. As we have seen, spectrally resolved measure- 
ments provide much better constraints on the actual transit 
depth than unresolved observations. We believe that it would 
be useful to repeat the O 1 transit depth measurements by us- 
ing this instrument. 

5. CONCLUSION 

Transmission spectroscopy in UV wavelength bands is a 
rich source of information about the upper atmospheres of ex- 
trasolar planets. In this paper we have introduced a generic 



method that can be used to interpret and analyze UV tran- 
sit light curves. This method is based on tracing the emit- 
ted flux from the stellar atmosphere through the atmosphere 
of the planet and the interstellar medium to the observing 
instrument at Earth. We have demonstrated the method in 
practice by applying it to the existing UV transit depth mea- 
surements of HD209458b in the H Lya and O 1 triplet lines 
(VM03,VM04,BJ10). In order to interpret the measured ti-an- 
sit depths we used a simple empirical model of the thermo- 
sphere based on the generic feat ures of more complex models 
(e.g.,|Y elle 2004; Garcia Mun oz|2007[ [Koskinen et al.|2010| l 
to simulate absorpFion by the occulting atmosphere. 

The H I transit depth is sensitive to interstellar absorption 
and the full-width transit depth reflects the optical depth of 
the atmosphere in the wings of the H Lya line. We found 
that it is easy to generate models that explain the observed 
absorption without the need to introduce external ENAs or 
hydrogen atoms accelerated by radiation pressure outside the 
Roche lobe. According to the best-fitting models (see models 
Ml and M7 in Table Hi, the mean temperature of the thermo- 
sphere is T = 8000-1X000 K and the H2/H dissociation front 
is located at pb = 0.1-1 /ibar {zb^ 1-1 Rp)- The upper bound- 
ary of the model is located at ^ 2.9 Rp, which is near the 
boundary of the Roche lobe. Below the density profiles 
are approximately hydrostatic and above the atmosphere is 
mostly ionized. 

By using the Ml model, we obtain a full- width H 1 transit 
depth of 6.6 % and the model transmission of the H Lya line 
matches the data points reported by BJ08. The apparent dis- 
agreement between BJ08 and the earlier analysis by VM03 
arises from differences in the treatment of the data and dif- 
ferent definitions of the transit depth. In particular, there is 
no definite observational evidence for a cometary tail follow- 
ing the planet in the transit light curve and the absorption is 
not significantly stronger in the blue side of the H Lya line. 
However, we agree with VM08 that the core of the absorption 
line is optically thick up to the Roche lobe and thus that the 
atmosphere is evaporating. We estimate a mass loss rate of 
dM/dt wlO^-lO*^ kg s"' based on the density profile of the 
Ml model and a range of possible radial velocities. 

It is possible that ENAs are present outside the Roche lobe 
as suggested by ElO, but the optical depth of the ENA clouds 
has to be consistent with the underlying models of the ther- 
mosphere. We do not agree with ElO that the optical depth 
of the thermosphere below z= 2.8 Rp is small. We do be- 
lieve, however, that one-dimensional hydrodynamic models 
are not adequate in modeling the distribution of ionized gases 
outside the Roche lobe. Multidimensional plasma models are 
required to study the interaction of the ionosphere with the 
stellar wind and magnetospheric plasma self-consistently. 

The mean O 1 transit depth of 10.5 % is only 1.93cr away 
from the FUV continuum transit depth of ^2 %. By us- 
ing the Ml model, we obtain an O 1 transit depth of 4.3 %, 
which is marginally consistent with the observations. Full- 
width O 1 transit depths that are comparable or slightly higher 
than the H 1 depth are possible if the abundance of oxygen 
in the thermosphere is supersolar, the atmosphere is escaping 
with supersonic velocities or if large external clouds of neu- 
tral oxygen are present above the Roche lobe. More precise 
measurements are required to constrain the O 1 transit depth 
further Out of the currently available instruments, the STIS 
G140M grating can be used to obtain repeated measurements 
with higher S/N and spectral resolution. 
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